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I give a short summary of jet definition algorithms and recent progress in the quanti- 
tative description of jet production. 
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. 1. Introduction 



Over the last decade detailed theoretical description of jet physics at high energy col- 
liders has been established at the next-to-leading order (NLO) accuracy [ i]- The 
^ I theoretical predictions have been successfully compared with the experimental data in 
the case of three jet production in e~^e~ annihilation, one or two jet production in photo- 
• production, deep inelastic scattering at HERA, and one or two jet production at proton- 



antiproton collisions at Tevatron 



D ! Jets are footprints of quarks and gluons produced at short distances. They are defined 
^ ' as some collections of final state hadrons such that the relative angular distances between 
the hadrons in momentum space are small. Important cancellation theorems valid in 
^ ■ all orders of perturbation theory suggest that infrared safe global jet observables can 
quantitatively be described in terms the weak perturbative dynamics of the point like 
partons of QCD. A jet observable is infrared safe if it does not change by adding or 
removing a soft particle from the jet or by splitting an ultra relativistic particle into two 
coUinear particles within the jet. 

The data with highest jet energy at Tevatron could give us the most stringent test of the 
QCD dynamics at short distances and allow for an efficient search for effects of deviations 
from the Standard Model predictions [p. The data of CDF and DO collaborations at 
Tevatron have about the same (< 20%) or better statistical accuracy as the theoretical 
predictions. At the highest jet energies, the uncertainties in the jet energy calibration and 
the errors of the fitted values of parton number densities, however, lead up to systematic 
errors of ~ 50%. 

In the case of three jet production in e~^e~, the experimental systematic errors are 
smaller and the theoretical uncertainties are considerably larger than the current experi- 
mental errors [ Therefore the theorists have strong motivations to further improve the 
accuracy of their predicitons. 

Recently, we could witness important theoretical progress in four areas of jet physics. 
First, the NLO calculations could be extended to 5-leg processes (four jet production in 
e~^e~ annihilation [ 0, |H, three jet production in hadron-hadron collisions [ ^, |l^ etc.). 
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Secondly, remarkable progress has been made towards the ambitious goal of calculating jet 
cross sections of 4-leg processes in next-to-next-to leading order (NNLO) accuracy [ pJ]] . 
In particular, the NNLO virtual corrections could be calculated analytically for four leg 
amplitudes. Thirdly, the NLO calculations could be improved in the threshold region by 
resumming the large logarithmic contributions to all order [ |12| . Finally new techniques 
have been developed to calculate many jet processes in leading order [ Below I discuss 
difficulties related to jet definitions, I briefiy review the general theoretical framework of 
the NLO calculations, I describe in some detail the new NLO results obtained for 5-leg 
processes, finally I summarize the main results on NNLO corrections. 



2. Jet definition algorithms 

Jet definitions are based on a selection and a recombination algorithm. The selection al- 
gorithm selects the group of particles which form the jet and the recombination algorithm 
specifies how to construct the kinematicl variables of the jet in terms of the momenta of 
the selected hadrons. 

The jet algorithms used in the data analysis are more complicated than the simple 
theoretical definitions in terms of one, two or three partons. Some 'auxiliary' features like 
the phenomena of merging and splitting, or the use of seeds are not always completely 
specified. These issues may infiuence the infrared sensitivity and the size of the hadronic 
fragmention corrections. In a careful analysis one has to require that i) the jet selection 
process, the jet kinematic variables, preclustering, merging, splitting, the role of the 
underlying evens are fully specified; ii) The fully specifed algorithm should be infrared 
safe, independent from the detector properties and the algorithms should be simple to 
use; iii) The algorithm has to be defined in the same way at parton, hadron and detector 
level (for a recent detailed discussion for hadron colliders see [ p^,?]). 

Although the algorithms have a large amount of arbitrariness they should fulfil a num- 
ber of important theoretical requirements, such as insensitivity to soft radiation and to 
coUinear splitting, or invariance under boost in the case of hadron colliders etc.. In the lat- 
ter case the use of transverse momentum and rapidities as kinematic variables is preferred. 
Another constraint is suggested by the study of resummation corrections. Resummations 
can only be carried out if the boundaries of the inclusive jet kinematic variables are de- 
fined independently from the number of final state particles [ |I6| . This prefers Lorenz 
covariant recombination schemes (such as the E-scheme). Finally, the algorithm should 
be as simple as possible. At hadron colliders we do not have one best algorithm. The 
cone algorithm is broadly used while the kj- algorithm is preferred more by the theorists. 

2.1. Cone jet algorithm for hadron-hadron collisions 

The cylindrical shape of the detectors and boost invariance suggest that the kinematics 
has to be described in terms rapidities and azimuthal angles [ 0, 0]. In the 2-dimensional 
1] X (f) lego-plot the hadrons or partons constituting the jets lie within a cone of radius R. 
The trial cone in the lego plot is given by its radius and the value of its center (77*-^ , 0^). 
The cone is adjusted in such a way that the geometric center of the cone agrees with the 
£^r-weighted recombined values of the particles within the cone. All particles within the 
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trial cone fulfil 

^^r^c_^;)2^^^c_^2)^R, (1) 

A stable cone (protojet) is obtained if the " physical" center of the cone defined by the 
recombination algorithm 

R _ Si ?7t-^T , R _ Si 0T-^T _ pi (r)\ 

Vt - — ^;r — ) 0T - — ^;r — ) - [^) 

T T i 

coincides with its geometrical center [j]^ ^cfP) = {ri^,(j)^). In this case it is natural to 
identify the jet variables with the recombined values of the stable cone {E^,rj'^ = 
{E^ , rj^ , (j)^) . After identification of the jet as group of particles within the stable cone 
one can construct the jet kinematical variables with using some recombination scheme, 
for example 

i 

E^ 

6^ = arctan^, 6'^ = arctan^, (4) 

7]^ = - In l^tan y ^ , E^ = E"^ sin 6^ (5) 

Different recombination algorithms give aproximately the same kinematical values if 
M'^ << E^. The boost invariant recombination scheme, however, is a better estima- 
tor of the jet variables. I have noted above that boost invariant variables are not suitable 
for resummation studies since the kinematic boundary of the Ej- depends on the number 
of partons. Resummation is consistent with e.g. the E-scheme which recombines the jet 
four momenta as = 

If the number of the final states particles are large the method of moving the trial 
cone gives a slow algorithm. One can speed it up with starting the cone iteration at 
the center of seed towers which passed a minimum cut-off. But the seeded algorith may 
have problems with sensitivity to the emission of soft gluons. The algorithm allows for 
jet overlaps, therefore, one should also specify the details of jet merging and jet splitting. 
This allows for many ad hoc options and it is difficult to ensure that the definition of jets 
are the same at the detector and the parton levels. In a quantitative analysis all these 
details should be treated with great care. 

2.2. algorithms 

The kx algorithm was designed to avoid jet overlaps and problems with kinematical 



boundaries in the case of resummation [ [T5|. I summarize the version of kx algorithm 



suggested by Ellis and Soper. The algorithm starts with the initial list of particles and 
an empty list of jets. Performing the algorithm we get an empty list of particles and the 
list of jets, each separated by . This regrouping of particles into jets goes iteratively 
through five steps. 

i) Every particle (pseudoparticle) i and every particle pair in the particle list is 
associated with a d— value 

~D2 



di=plT dij = minipl^, plj.)—f- (6) 
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where D is a. free parameter (the usual choice for its value is D = 1) and Ai?^^ is 
the square of the distance between the particles in the lego-plot '^Rfj = — rjjY + 

ii) Find (imin = min((ij, dij). 

iii) If (ijnin = dij replace the particle pair {i^j) with a pseudoparticle. Its momentum 
Pij is calculated via the rules of the recombination scheme. For example using the 
i?-scheme pij = pi +Pj- 

iv) If (imin = di remove particle (pseudoparticle) i from the list of particles and add it 
to the list of jets. 

v) If at this step the particle (pseudoparticle) list ist not empty got to step i). 

In the experimental analysis the algorithm has to combined with some preclustering 
procedure in order to keep the jet analysis computionally feasible and to diminish the 
detector dependence of the alghorithm. The algorithm has the tendency to reconstruct 
more energy from calorimeter noise, pile-up and underlying event and multi pp interactions 



than a cone algorithm [ |14[. The jet momentum resolution appears, however, to be the 
same for kx jets and cone jets. Therefore, it appears that the use of kjn jets has clear 
overall advantages. 



3. NLO cross sections 
3.1. Formalism 

The lowest order cross-sections strongly depend on the unphysical renormalization and 
factorization scales. Higher order cross sections reduce this sensitivity with a factor of 
0{as)- The calculation of the NLO corrections is technically involved since the virtual 
and gluon bremsstrahlung contributions are separately divergent ( soft and collinear sin- 
gularities). Fortunately, the divergent pieces are universal in the sense that they are 
given by the Born cross sections times a universal process independent singular factor [ 
H, ^. The singular parts of the virtual cross sections can be cancelled analytically with 
contributions of of the real contributions with the use of local counter terms. They are 
used to subtract the real contributions in the singular soft and collinear regions. This 
subtraction procedure is consistent with numerical Monte Carlo evaluation of the phase 
space integrals provided one calculates physical quantities which are defined in terms of 
infrared safe measurement functions specified below. In the case of the collisions of hadron 
A with hadron B, the physical cross section is given in terms of parton number densities 
fh/nix,^) and parton-parton scattering cross sections CTa^b as follows 

(^{PA,Pb) = Ea,b I dXadXbfa/AiXb, IJ-F)fb/BiXb, fip) (7) 
X {^a!b{^aPA,XbPB) + O'^fi'^ {XaPA, XbPs)] 

where for example the LO cross section for the physical quantity S is given by the phase 
space integral over the squared matrix element of the 2-to-n process, weighted by the 
appropriate measurement function S{pa,Pb',Pi, ■■■Pn) 

a^a?{pa,Pb-AS^n)]) = Jja^{Pa,Pb) = J dT^''^\Ma,b\'S^'^\pa,Pb,PU..Pn) 
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Inclusive Jet cross section 
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Figure 1. Inclusive jet cross-section [ ^ as 
function of the jet energy at various pseudo 
rapidity intervals. The CDF data are com- 
pared with the predictions calculated by the 
JETRAD NLO MC program [m. 



Figure 2. Inclusive jet cross- 
section as function of the trans- 
verse energy measured by the 
CDF Collaboration [0. The 
solid curve is the NLO QCD pre- 
diction \ 131. 



The finite parton cross sections in NLO are obatained by summing the virtual, real and 
counter term contributions 

^a'^^iPa^Pb) = I ^a,biPa^Pb) + [ ^lb{Pa,Pb) + / ^aliPa.Pb) 

The individual contributions are evaluated in d-dimension. In the numerical evaluation 
of the phase space integral over the real part we need to subtract the singular contribu- 
tions locally. Fortunately there are several general methods for constructing such local 



subtarction terms analytically [|^, y, |T6[. We can write 

= {Kb - Kl} + Kb + Kb + Kl ' (8) 

The local subtraction term da^°^ is a suitable approximation of da^^,. They become equal 
in the singular regions. In order to cancel analytically the singular parts of the virtual 
contributions we must be able to carry out analytically the integrals over the single par- 
ton subspaces in ci = 2 — 2e dimension. After performing this integration over the local 
subtraction terms the singular pieces of the terms in eq. can be cancelled analytically 
and the remaning part can be evaluated numerically in four dimensions. The cancellation 
mechanism of the soft and coUinear singularities of the virtual corrections against the sin- 
gular part of the real contribution is independent of the form of the measurement functions 
provided that they are insensitive to coUinear splitting and soft emission. That is in the 
soft or/and collinear configurations the measurement functions must fulfil the condition 

5("+^Hpa,..Pn+l)=5("nPa,...p„) (9) 
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Figure 3. Comparison of the data and the 
ory in the case of inclusive jet production, 
are defined by the the kr algorithm [ |17| . 



Figure 4. Theoretical prediction for 
inclusive three jet production at the 
Tevatron ( Z. Nagy [0). The band in- 
dicate the theoretical uncertainty due 
to the variation of the renormalization 
and factorization scales xr^f between 
0.5 and 2. The inset shows the K- 
factor. 



The existence of universal local subtraction terms [ ^, [T6| is crucial for the method. 
The most widely used implementation is the one by Catani and Seymour [ |16|. The 
construction of the NLO parton level Monte Carlo programs for 4-leg processes is by now 
a routine (although rather laborious) work. The success of the NLO description of jet 
production at hadron colliders is well illustrated by the plots shown figures 1 and 2. We 
see a spectacular agreement in a wide range of kinematical variables. The ambiguities 
due to errors in the parton number densities are shown in figure 3. 



4. NLO description of 5-leg process 

4.1. New results for e+e^ 

The final experimental analysis of 4-jet production and its comparison with the NLO 
theoretical predictions was recently completed [ 0. The most remarkable result is that 
we obtain a more precise measurement of the strong coupling constant from the analysis 
of the 4-jet data than from the 3-jet data. This is not completely suprising since the 
cross section of four jet production is more sensitive to than the three jet cross section, 
furthemore due to the large number of events the statistical accuracy is not the dominating 
the experimental error. 
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Figure 5. Comparison of the scale depen- 
dence of the three jet cross sections ob- 
tained in LO and NLO accuracy [ 
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4.2. New results for pp{p) and ep scattering 

Important progress in this topics is that the theoretical results [ ^ (where the singular- 
ities are not manifestly cancelled) could be implemented in a efficient parton level C++ 
Monte Carlo program (called NLO-I--I- ) both for ep [ and for pp{p) [ ^ scattering. 
Similarly to the case of e^e~ annihilation, the accuracy of the theoretical prediction is 
greatly improved. This is illustrated in figure 4 and 5. In fig. 4 inclusive jet cross section 
for three jet production is plotted using the algorithm. As we mentioned above, in 
multi jet production the use of the cone algorithm is cumbersome. It is assumed that the 
jets are produced in the pseudo rapidity interval |r/| < 4 and at energy of y/s = 1800 GeV, 
with minimum transverse jet energy Et > 50 GeV. The cross section is plotted as a func- 
tion of the transverse energy of the leading jet The hard scattering factorization 
and renormalization scale is chosen to be The theoretical ambiguity is indicated as 
a band when the hard scattering scale is changed by a factor in the range 0.5 < xfi,f < 2. 
One can see that at NLO the cross section values show much less sensitivity to the value 
of the hard scattering scale as the LO cross sections. Hopefully, these predictions can be 
compared with the Tevatron data in the near future. In fig. 5 integrated cross section 
values of three jet production (with the same kinematics as in fig. 4 but integrated over 
E^j^^ > 100 GeV) are plotted function of the ratio jJ^/Elj^^ where n is the value of the 
hard scattering scale. 

Due to crossing symmetry the e~^e~ NLO four jet matrix elements allow also the evalu- 
ation of the DIS 3+1 jet processes at NLO. The implementation of the matrix elements [ 
into C++ Monte Carlo program with the corresponding local subtraction terms has 
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been recently carried out [ jTO) . The new result is crucial in the quantitative comparison of 
the data with the theory. As one can see in fig. 6, at NLO, we get substantial reduction in 
the precision of the theoretical prediction and at the same time the large NLO corrections 
are needed to bring the data in agreement with the theory. 



5. NNLO calculations 



In the last few years, the most spectacular progress has been achieved in developing the 
technics of calculating cross section values at NNLO accuracy. This requires the analytic 
computation of two loop corrections to 4-leg QCD processes (for a recent more complete 
review see The basic ingredient behind this success is the use of two technical tricks. 

First, if one uses integration by parts [|2^ and Lorenz invariance [ 21] identities one can 
reduce the very large number of different integrals into a few master integrals. This step 
can be very well computerized and by now several computer algorithms are available which 
perform automatically this reduction. Second, all the remaining master integrals (double 
box diagrams, etc.) could be evaluated analytically. Three different methods have been 
used: i) Mellin-Barnes integral transformations [ ^ ; ii) inhomogeneous linear differential 

p3| . In addition. 



equations [ ^ and iii) expansions in terms of nested harmonic sums 
integrals could be calculated numerically using sector decomposition 



, the 
Very recently, 

m 



New 



the differential equation method could be extended also to phase space integrals 
classes of functions appear in the final solutions. Numerical algorithms are available for 
the fast evaluation of these functions. All parton-parton scattering amplitudes could be 
calculated analytically at two loop order. Analytic expressions of the virtual corrections 
to the cross sections of three jet production in e~^e~ annihilation have been published 
recently [ As a first phenomenological application Bern et.al. [ ^] have computed 
the cross section of two photon production at LHC in NNLO accuracy. In the case of this 
process the treatment of the bremsstrahlung corrections are considerable easier than in 
the general case. For processes involving initial hadrons, the parton evolution has to be 
evaluated also at NNLO accuracy. The results for the non-singlet case have been published 
recently by Moch and Vermaseren. This is an outstanding technical achievement. If the 
NNLO splitting functions 
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will be known the fitting of the initial parton density 
functions can be carried out at NNLO accuracy. Depending on the nature of the Beyond 
the Standard Model Physics, the newly achieved theoretical results may be crucial for 
the success of the experimental programs at future collider experiments. The impact of 
these developments can not be overestimated for the future high precision measurements 
at HERA, Tevatron, LHC and TESLA. 
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